Vacuum Rabi oscillation induced by virtual photons in the ultrastrong coupling regime 
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We present an interaction scheme that exhibits a dynamical consequence of virtual photons carried 
by a vacuum-field dressed two-level atom in the ultrastrong coupling regime. We show that, with 
the aid of an external driving field, virtual photons provide a transition matrix element that enables 
the atom to evolve coherently and reversibly to an auxiliary level accompanied by the emission of a 
real photon. The process corresponds to a type of vacuum Rabi oscillation, and we show that the 
efi^ective vacuum Rabi frequency is proportional to the amplitude of a single virtual photon in the 
ground state. Therefore the interaction scheme could serve as a probe of ground state structures in 
the ultrastrong coupling regime. 

PACS numbers: 42.50.Pq, 42.50.Ct, 42.50.Lc 



A single-mode electromagnetic field interacting with a 
two-level atom has been a fundamental model in quan- 
tum optics capturing the physics of resonant light-matter 
interaction. In particular, the Jaynes-Cummings (JC) 
model P, [2I , which describes the regime where the inter- 
action energy HX is much smaller than the energy scale 
of an atom Huja and a photon HlUc, has tremendous ap- 
plications in cavity QED |3| and trapped ion systems 
Q . Recently, there has been considerable research inter- 
est in the ultrastrong coupling regime where A becomes 
comparable to ujc and uja- Such a regime has been ex- 
plored by experiments in various related systems with 
artificial atoms and cavity photon resonators, including 
superconducting qubit in coplanar waveguide [6| or LC 
resonator JJ , microcavities embedding doped quantum 
wells 0, Q, and two-dimensional electron gas coupled 
to metamaterial resonators In addition, theoreti- 

cal investigations have also found novel phenomena in 
the ultrastrong coupling regime, such as the asymmetry 
of vacuum Rabi-splitting [ill, photon blockade [l^] , non- 
classical states generation [ij], superradiance transition 
[3], and collapse and revivals dynamics flst . 

A key feature in the ultrastrong coupling regime is the 
significant number of virtual photons existing around the 
vacuum-field dressed atom. These virtual photons are 
generated by counter-rotating terms in the Hamiltonian, 
and they can have direct physical consequences. For ex- 
ample, by modulating the atom-field coupling strength 
virtual photons can be released as a form of quantum 
vacuum radiation [l6| . In this paper we address a dif- 
ferent effect of the vacuum-field dressed atom, namely, a 
kind of vacuum Rabi oscillations that would not occur if 
virtual photons are absent. 

Specifically, we investigate the quantum dynamics of 
a driven quantum Rabi model. The configuration of our 
system is shown in Fig. 1 in which a S-type three-level 
atom is confined in a single-mode cavity. The atomic lev- 
els \g) and |e) are coupled to a cavity field of frequency 
ujc- These two atomic levels and the cavity field mode 
constitute a Rabi model. In addition, there is an ex- 
ternal classical field driving the transition between |e) 
and the third atomic level |/). We note that some the- 



FIG. 1: (Color online) Interaction scheme of a H— type three- 
level atom in a cavity. The atomic states \g) and \e) and a 
cavity field mode of frequency ujc form a quantum Rabi model 
described by Hr, and an external classical field of frequency 
Ljp drives the transition between |e) and |/). 



oretical aspects of three-level artificial atoms in circuit 
QED was discussed in [17[, and S— type superconducting 
atoms have been demonstrated in experiments (l8l - [20| . 
Recently Carusotto et al. have studied the dynamics of 



a related system in a different driving configuration [21| . 
The Hamiltonian of our system is given by (h = 1), 

H = HR + iOf\f){f\ + ncos^pt{\f){e\ + \e){f\) (1) 

where Hn is the Hamiltonian of the Rabi model [221, 



H 



^0 , 



{\e){e\^\g){g\)+u,a^a+\{a+a^){\g){e\ + \e){g\). 

(2) 

Here is the (bare) transition frequency between |e) and 
\g), and w/— wo/2 is the transition frequency between |/) 
and |e). The parameter A denotes the atom-cavity cou- 
pling strength, and the classical driving field has a fre- 
quency Up and an interaction strength il. In writing HR^ 
we have kept counter-rotating terms because A is com- 
parable to iOc in the ultrastrong coupling regime. Note 
that the coupling between the cavity mode and the level 
I/) is assumed to be weak and so that it is not included 
in the Hamiltonian. 

Initially the system is prepared in the ground state of 
Hn, which is the lowest-energy state of the system in the 
absence of the driving field. Our task is to determine the 
dynamics after the driving field is turned on. To analyze 
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the problem, we apply a unitary transformation to sim- 
plify the Hamiltonian. It is known that for low energy 
states of the Rabi model, Hn can be transformed to into 
a form of Jaynes-Cummings Hamiltonian approximately 
by a unitary operator e~'^ ^231]. Here the operator S and 
its parameters are defined by: 

5=-(l5)(e| + |e)(5|)(at-a), (3) 



where f2' = •q^/'^^ is a renormalized driving field strength, 
and the last term indicates a new coupling between \g) 
and I/) through the cavity field mode. 

A further simplification can be made by exploiting res- 
onance when Up is tuned to a certain resonance frequency 
defined by the undriven system. In this paper we con- 
sider the resonance at 



•q — exp( 



)• 



(4) 
(5) 



(11) 



Then it can be shown that H'j^ — e^Hjf e ^ is approxi- 
mately given by [H-iil 



UJ, 



r(|e)(e| - \g){9\) + ^c«ta + A'(a|e)(5| + a^\g){e\) 



+ ^{^-2){\e){e\ + \g){g\) 

UJc 

Hjc 



which corresponds to the transition between \g, 0) to 
/, 1), since the square bracket term is the approxi- 
mate ground state energy of Hji by the transformation 
method. By the condition (11), \g,0) and |/, 1) are res- 
onantly coupled, but |/, 1) and |e, 1) is far away from 
resonance (the corresponding detuning is of order ujc). 
Therefore if J7' is not too strong, the system is confined 
to the two resonantly coupled states, i.e., all off-resonant 
transitions may be ignored. In this way H' in the inter- 
action picture is reduced to 



(6) 



where Hjc describes a JC model in which the atomic fre- 
quency and cavity-atom interaction strength are renor- 
malized as uj'q = rywo and A' — 2r]UJo^X/ujc, respectively. 

Note that Hjc in Eq. (6) is an approximation to H'j^, 
and the difference — Hjc describes multi-photon pro- 
cesses that correspond to higher order corrections |23l - 
Since \g, 0) is the ground state oiHjc, e~^\g, 0) is an 
approximated ground state of Hn in the original frame. 
The accuracy of such an approximation has been tested 
in Ref. [l^. Specifically, if A is comparable but smaller 



l/,l)(5,0|). (12) 



Eq. (12) indicates that the system would execute a form 
of vacuum Rabi oscillations, in which |g, 0) behaves as an 
excited atom in the vacuum field, and |/, 1) behaves as 
an ground atom with a single photon. In cavity QED, 
such oscillations lead to vacuum Rabi splitting j27rt29i] . 
Note that the effective vacuum Rabi frequency here is 
A^il'/wc, which is significant in the ultrastrong coupling 
regime where A is comparable to uJc- 

It is useful to go back to the original frame in which 



than UJc, the ground state energy of Hjc has a good ^he Rabi oscillations occur between the states e-^|/,l) 
agreement with that of Hr obtained by exact numerical and e-^|g, 0). Since e"^]/, 1) = |/, 1), an initial ground 
calculations over a range of parameters. For example in 
the case Wc = wq = 2A, the approximated ground state 
energy obtained by Hjc has the percentage error about 
0.65%. 

Now we perform the transformation for our system 
Hamiltonian i7, which becomes. 



H' = e^He-'^ 

« H.jc+ujf\f){f\ 

+^lcosujpt{e^\e){f\ + \f){e\ 



(7) 

Since e^|e) = cosh[^(at - a)] |e) -|-sinh[^(a^ - a)]|g), we 
expand the hyperbolic sine and cosine operator functions 
in normal order up to first order in A^/wq, 



cosh 
sinh 



(at - a) 



^1/4 



(8) 



^(at-a) 

UJc 



r^m^^^a^-a) (9) 

Uc 



Therefore the transformed Hamiltonian becomes, 
H' « Hjc+ujf\f){f\ + n'cosujpti\f){e\ + \e){f\) 

+ ^n'cosu;pt{\g){f\ - \f){g\) [a^ - a) (10) 



state will evolve to |/, 1) after half of a Rabi period. If 
we switch off the external field at this moment, the single 
photon described by |/, 1) will be free to escape the cavity 
because the atom in the state j/) does not couple to 
the cavity field when 51 = 0, i.e., the photon cannot be 
reabsorbed by the atom. In this way, a tt pulse of the 
driving field can generate a real photon deterministically 
while the atom is excited to the |/} state. 

To gain a better insight of the physical process without 
relying on the approximation made in Eqs. (6) and (10), 
we express the Hamiltonian by the eigenbasis of Hr. Let 
\ipn) be an eigenvector of Hr with the eigenvalue A„, i.e., 
HrI^ii) = Xnl'ipn) (the ground state is denoted by IV'o)), 
and consider the expansion \e,n) — J2m'^nm\'4'm) with 
the coefficients Cnm = (V'mlGi"')- Therefore 

l/)(e| ^Y(^nm\f,n){llJm\- (13) 

n nm 

In this way, the Hamiltonian (1) in the interaction picture 
becomes. 

Hi = f7cosc.pt^e'('^^+""^-^-)*<„J/,n)(V'™K/i.c. 

nm 

(14) 
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FIG. 2: (Color online) Probability amplitude of je, 1) in the 
ground state of Hr as a function of the coupling strength A 
for the uq = uj^ case. The solid red line corresponds to exact 
numerical values, and the dashed blue line is obtained from 
the approximated ground state e~'^\g, 0) according to Eq. (6). 

At the resonant frequency ujp — ujf + ujc — ^Oi IV'o) a-nd 
I/, 1) are resonantly coupled. If we keep only the resonant 
terms, then we have 

Hj^^\f,l){^o\+h.c. (15) 

Comparing with Hj in Eq. (12) and noting that |i/'o) ~ 
e~^\g,0), Hi describes the same type of resonant inter- 
action as H'j. However, we emphasize that Hj in Eq. 
(15) is a more accurate interaction Hamiltonian than H'j 
because Hj is derived directly from the eigenbasis of Hj^ 
without making use of the approximation in Eq. (6). In 
this sense, the resonant condition (11) can be improved 
by replacing the square bracket term by Aq. 

The role of virtual photons is now explicitly seen in Eq. 
(15) through the effective vacuum Rabi frequency f2|cio|. 
This is because cio is precisely the probability amplitude 
of a single virtual photon state in ji/jo). In other words, 
we may interpret that the interaction described in Eq. 
(15) is induced or mediated by a virtual photon. In Fig. 
2, we plot cio (solid line) as a function of X/ujc for the 
case ujc — ujQ, and the figure shows that the magnitude of 
cio is appreciable in the ultrastrong coupling regime. As 
a comparison, we also plot the approximate amplitude 
Cio ~ —Tj^^^^.X/^c (dashed line) obtained from e~"^|g, 0). 
For the parameters used in Fig. 2, we see that the approx- 
imation agrees well with the exact numerical calculation 
up to X/ojc < 0.6. 

We have tested our prediction of the virtual-photon- 
induced Rabi oscillations by solving numerically the 
Schrodinger equation defined by the Hamiltonian (1) 
with the initial state {ipo). In Fig. 3 we plot the exact nu- 
merical probability Pi / of the system in the state |/, 1) as 
a function of time. The parameter A = cjc/2 used in the 
figure is served as an example of ultrastrong coupling. We 
see the Rabi cycles as predicted by the Hamiltonians (12) 
or (15) for relatively weak driving fields with ft < OAlUc- 
At a stronger driving field with fi = O.Swc (red solid line), 
and there is a high frequency pattern due to counter ro- 
tating terms of the classical driving field, and the Rabi 
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FIG. 3: (Color online) Probability of |/, 1) as a function 
of time for fl = 0.2a;c(blue long dashed), 0.4ljc (green short 
dashed) and O.StOc (red solid). The parameters used are: 
A — 0.5oJc, oJc = i^o ~ '^//S, ojp = ujf + iOc — Aq, and the 
numerical ground state energy Ao — — 0.633tJc. The figure is 
essentially the same if uup in Eq. (11) is used. 



oscillations are less perfect in the sense that the maxi- 
mum Pij K 0.9 is smaller than one. Such a behavior 
is understood because the off-resonance transitions ne- 
glected in Eq. (12) or (15) would generate energy shifts 
which in turn could bring the driven system out of reso- 
nance. As a result, the amplitude of oscillations in Pi/ is 
reduced. Since these energy shifts are generally propor- 
tional to Sl^ , as long as O is small compared with detun- 
ings associated with off-resonance transitions, it would 
be safe to use Eq. (15), and this is demonstrated in Fig. 
3 for up to QAoJc- 

Finally, it is worth noting that the Hamiltonian in Eq. 
(14) has higher resonances ai ujp — u f + nujc — Aq for odd 
positive integers n. The requirement of an odd n is be- 
cause iV'o) has a definite parity in which the atomic state 
|e) and odd photon numbers are connected. In the case 
n = 3, the driving field at the corresponding uip would 
resonantly excite the atom to |/) with the emission of 
three real photons. The effective Hamiltonian would be 
of the same form of (15), but with |/, 1) and c*g replaced 
by I/, 3) and Cgg, i.e., the effective Rabi frequency is pro- 
portional to |c3o|. Such a three-photon resonance was 
also observed in our numerical calculations. 

To conclude, we have shown that virtual photons in 
the ultrastrong coupling regime can play a key role in 
quantum dynamics by providing the transition matrix el- 
ements that allow the system to access relevant quantum 
states of interest. In our scheme, the system can exhibit a 
form of vacuum Rabi oscillations which can be considered 
as a signature of virtual photons. Since our main focus in 
this paper is on the interaction induced by virtual pho- 
tons, decoherence effects have not been included in the 
discussion. However, as long as the decoherence times 
is sufficiently short, coherent dynamics predicted by the 
Hamiltonian (12) or (15) would be justified. Specifically, 
given a vacuum Rabi period T « 2tt lo c / X^Ct' , the cav- 
ity field damping rate 7c and atomic decay rate ja, the 
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condition ^jT <^ 1 {j = c,A) ensures that the system 
can execute a Rabi cycle without being afTected by the 
damping, and this is achievable in the ultrastrong cou- 
pling regime with moderate small 7's. For the parameters 
used in Fig. 3, for example, 7^ < 10~^cJc would be suffi- 
cient. We emphasize that a finite interaction time within 
T is of practical importance, since the interaction (f2) or 
(f 5) is switchable via the driving field. This feature could 
be a tool for performing quantum operations on qubits 
formed by the atom or the field, as well as for determin- 
istic single-photon generation j30l - l33 |. In addition, since 
the effective vacuum Rabi frequency is proportional to 



the corresponding virtual photon amplitude, our scheme 
can be used to probe the ground state structure of the 
quantum Rabi model. 
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